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Finite temperature transition in lattice QCD with degenerate Wilson quarks is investigated on an Ni = 4 
lattice, using a renormalization group improved gauge action. We find the following for the Np = 2 case: 1) The 
transition is smooth for a wide range of the quark mass. 2) The chiral transition is continuous. 3) The chiral 
condensation well satisfies a scaling relation with the critical exponents of the 3 dimensional 0(4) spin model. 
For Np = 3, we find that the chiral transition is of first order. 



1. Introduction 

In our previous work [1] we investigated fi- 
nite temperature lattice QCD with two degener- 
ate Wilson quarks using a renormalization group 
improved gauge action [2] and the standard Wil- 
son quark action. We showed that the unex- 
pected sharp transition observed with the stan- 
dard gauge action in the intermediate mass range 
[3] disappears on an TVt = 4 lattice with the im- 
proved action in the quark mass region where the 
sharp transition was observed. In this article we 
extend this work to a wider range of /3 and ex- 
tensively investigate the dependence of various 
physical quantities on {3 and quark mass. Among 
others we check a scaling relation suggested by a 
universality argument [4]. We also investigate the 
three flavor case with the improved action. 

Numerical simulations for Np = 2 are done on 
8^ X 4 and 8* lattices. The critical hopping pa- 
rameter Kc is determined by rre^ = on the 8^ 
lattice and the lattice spacing a from mp at Kc- 
For Np = 3 we simulate 8^ X 4 and 12^ X 4 lattices. 

2. Nf=2 

In Fig. 1-a we show the Polyakov loop obtained 
with our improved action for TYj = 4 as a function 
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of l/K — l/Kc (Kc is interpolated when data is 
not available). The change of the Polyakov loop 
at the finite temperature transition/crossover Kf 
is very smooth for a wide range of f3. As we de- 
crease /3, Ki moves toward Kc monotonously and 
the change of the Polyakov loop becomes steeper. 
These properties are in sharp contrast with those 
observed for the standard action (Fig. 1-b). 

Fig. 2-a shows the pion screening mass squared 
TO^. As the quark mass decreases, devi- 
ates gradually from that in the low temperature 
phase. Cusps observed with the standard ac- 
tion[3] (Fig. 2-b) disappear for the whole range 
of (3. We also observe that the envelope of rre^ on 
the Nt = 4 lattice agrees excellently with for 
Nt=8. 

It is known that the quark screening mass 
defined by an axial-vector Ward identity [5.6] ex- 
hibits peculiar behavior with the standard action 
in the high temperature phase at /3 < 5.3: for 
Nf = 4 shows a sharp bend at Kf and the value in 
the high temperature phase does not agree with 
that in the low temperature phase[3,l]. As shown 
in Fig. 3. rUq with the improved action shows no 
such strange behavior: for N-t = 4 changes 
smoothly from the low temperature phase to the 
high temperature phase and agrees well with that 
obtained on the N-t = 8 lattice. 

It was pointed out[7] that the unexpected sharp 
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Figure 1. Polyakov loop for Nf = 4 versus 1/K — 
1/Kc obtained witli (a) the improved action and 
(b) tlie standard action. /3 = 2.1 (1.6) for the 
improved action corresponds in terms of a^^ to 
P = 5.1 (4.5) for the standard action (a~^ ~ 1.2 
(0.83) GeV). 



transition in the intermediate quark mass region 
with the standard action is probably due to the 
unusual relation between the line and the Kf 
line: the Kt line which deviates from the line 
at the chiral transition point approaches toward 
the Kc line again at /3 ^ 5.0 where the sharp 
transition was observed. In the case of the im- 
proved action, as shown in Fig. 4, the distance 
between Kc and Kt decreases monotonously to- 
ward the crossing point, i.e. the chiral transition 
point (3ct- In accord with this, the sharp tran- 
sition disappears with the improved action, as 
shown above. Thus we conclude that the unex- 
pected sharp transition observed with the stan- 



Figure 2. The same as Fig. 1 for rre^. 



dard action is a lattice artifact. 

Applying the method described in R.ef. [7.8], 
we have performed simulations on the line in 
order to locate fid and to determine the order 
of the transition there. In fig. 5 we show ?r),^ on 
the Kc line together with that on the Kt line as 
functions of p. The rre^ decreases toward zero 
as f3 decreases and no two state signals are ob- 
served. This implies that the chiral transition is 
continuous. A linear extrapolation of rre^ gives 
Pet ~ 1-4. The smoothness of physical cpiantities 
on Kt at intermediate masses shown in Figs. 1-a 
and 2-a strongly suggests that the transition is a 
crossover there. This is in accordance with the 
theoretical expectation that a continuous chiral 
transition turns into a crossover at finite m^. 

A universality argument suggests that, if the 
chiral transition in Np = 2 QCD is of second or- 
der, then it belongs to the same universality class 
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Figure 3. The same as Fig. 1-a for 2mq. 



as the 3 dimensional 0(4) Heisenberg model [4,9]. 
We can use the resulting scaling relations to make 
clear the nature of the chiral transition. Smooth 
behavior of physical quantities with the improved 
action, in sharp contrast with the case of the stan- 
dard action, encourages us to test the scaling with 
Wilson quarks. Taking correspondences of /3—/3ct, 
2rriq and the chiral condensation ('^V') in QCD to 
the reduced temperature t, the magnetic field h 
and the magnetization M, we test the scaling re- 
lation [4] 

where / is a scaling function and the /3 and the 6 
are critical exponents. (This scaling relation has 
been tested for staggered quarks using different 
correspondences of scaling variables [10].) 

The explicit chiral violation of Wilson fermions 
leads us to use a subtracted chiral condensation 
{■ipip) defined by a chiral Ward identity [5]: 

X 

where Psix) is the pseudo scalar density. (We 
have dropped the overall renormalization coeffi- 
cient which plays no role in our study of scaling 
properties.) We vary P^t as a parameter to find 
a value where the scaling relation is best satis- 
fied with the 0(4) critical exponents [11]. Us- 
ing all data in the range oi (3 = 1.6 - 2.0 and 
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Figure 4. Phase diagram for Np = 2 and Np = 3 
with the improved action. 

< 2mq < 0.9, we find that data are nicely on a 
universal curve when we set /3ct = 1.34(3). Fig. 6 
shows the scaling function Mhr^/^ as a function 
of the rescaled temperature th~^^^^. We note 
that this value of (irj, is consistent with our results 
for 77),^ shown in Fig. 5. because the universality 
predicts that on the both lines should bend 
upward when /3 approaches toward /3ct [4]. On the 
other hand, when we fix the exponents to mean 
field values, the data for {'4''>p) are not on a univer- 
sal curve so well as in the case of 0(4) exponents. 
Furthermore the best fit value [ict = 1.48(3) in 
this case is apparently too large when we consider 
the results in Fig. 5. Therefore we conclude that 
our data are consistent with the picture that the 
chiral transition belongs to the universality class 
of the 3-dimensional 0(4) Heisenberg model. The 
success of our scaling test strongly suggests that 
the chiral transition is of second order for Np = 2 
in the continuum limit. The direct extraction of 
each critical exponent and the comparison with 
that of 0(4) model are necessary to make this 
conclusion more definite. 

3. Np=3 

The phase diagram obtained so far has been 
given in Fig. 4. Simulations are made on the 
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Figure 5. on the Kc line and that on the Kt 
line versus (3 for Np = 2. 



line at /3 = 1.1,1.2 and 1.3 on a 12^ X 4 lat- 
tice. Using the method described in Ref. [7,8], 
we find 1.1 and observe two-state signals at 

P = 1.1. Therefore we conclude that the chiral 
transition is of first order. Simulations for finite 
quark mass are in progress. 
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Figure 6. A test of a scaling relation for the chiral 
condensation. See text for details. 
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